AD-A070  599 


UNCLASSIFIED 


MARYLAND  UNIV  COLLEGE  PARK  DEPT  OF  PHYSICS  AND  ASTRONOMY  F/G  20/5 
NONLINEAR  DEVELOPMENT  OF  THE  NEGATIVE-MASS  INSTABILITY  AND  ITS  — ETCCU) 
1975  R C DAVIDSON*  Y GOREN.  H UHM  N00014-75-C-0309 


8 79 


iN2A0?0599 


NONLINEAR  DEVELOPMENT  OF  THE  NEGATIVE -MASS  INSTABILITY  AND  ITS 


EFFECT  ON  INTENSE  MICROWAVE  GENERATION 

Y.  Goren  and  H.  Uhm 
Department  of  Physics  and  Astronomy 
University  of  Maryland,  College  Park,  Maryland  20742 

R.  C.  Davidson^ 

Division  of  Magnetic  Fusion  Energy 
Department  of  Energy,  Washington,  D.  C.  20545 

ABSTRACT 

The  nonlinear  development  of  the  negative-mass  instability 
is  investigated  for  a thin  E-layer  (A<<Rq)  located  inside  a 
cylindrical  waveguide.  A quasilinear  analysis  shows  that 
the  instability  saturates  by  the  development  of  an  energy 
spread  and  a corresponding  increase  in  thickness  A of  the  E- 
layer.  The  efficiency  of  microwave  generation  by  this 
instability  is  also  calculated. 
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1.  INTRODUCTION  AND  ASSUMPTIONS 


The  negative-mass  Instability  has  received  considerable 
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attention  as  a mechanism  for  intense  microwave  generation. 

For  the  most  part,  however,  analyses  of  this  instability  have  been 

limited  to  the  linear  regime.  The  purpose  of  this  paper  is  to 

investigate  the  nonlinear  development  of  the  negative-mass  instability, 

with  particular  emphasis  on  the  implications  for  intense  microwave 

generation.  The  analysis  is  carried  out  for  a thin  (A<<Rq) , 

relativistic  nonneutral  E-layer  aligned  parallel  to  a uniform  axial 

magnetic  field  B„e  and  located  inside  a grounded  cylindrical 
trV.2 

waveguide.  The  nonlinear  consequences  of  the  negative-mass  instability 
are  investigated  within  the  framework  of  a simple  quasilinear  model 
that  makes  use  of  total  energy  balance  for  the  system.  One  of  the 
most  important  features  of  the  analysis  is  that  the  instability 
saturates  by  the  development  of  an  energy  spread  and  a corresponding 
Increase  in  thickness  A of  the  E-layer.  The  corresponding  efficiency 
of  microwave  generation  is  also  calculated. 

As  illustrated  in  Fig.  1,  the  present  analysis  is  carried  out  for 
an  unneutralized  E-layer  aligned  parallel  to  a uniform  external  magnetic 
field  and  located  inside  a cylindrical  conducting  waveguide 

with  radius  Rc.  Cylindrical  polar  coordinates  (r,0,z)  are 
introduced,  and  the  main  assumptions  are  enumerated  below. 

(a)  The  E-layer  is  thin,  i.e.,  A/R0«l,  where  A-Rj-Rj^  is  the 
radial  thickness  and  Rq=*(R^+R2) /2  is  the  mean  radius  of  the  E-layer. 

(b)  The  electron  motion  is  ultrarelativistic  (yq>>1)  and  the 

mean  equilibrium  motion  of  the  E-layer  is  in  the  azimuthal  g 

direction  (i.e.,  V®-0,  where  V®  is  the  mean  axial  velocity  of  an 
z z 


electron  fluid  element). 
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(c)  It  is  further  assumed  that  v/yq  «1,  where  v=Nee  /me 

is  Budker's  parameter,  Ng  is  the  number  of  electrons  per  unit  axial 

2 

length  of  the  E-layer,  c is  the  speed  of  light  in  vacuo,  YQmc  is 
the  azimuthal  electron  energy  at  equilibrium  radius  Rq,  and  -e  and 
m are  the  charge  and  rest  mass,  respectively,  of  the  electron. 


2.  theoretical  model 


The  equilibrium  and  negative-mass  stability  properties  have  been 

investigated  in  Refs.  7-9  for  the  choice  of  a sharp-boundary  equilibrium 

in  which  the  electrons  h-ve  a rectangular  density  profile,  i.e.,  n^(r)« 

N / (2ttR.A)  =const . , for  R,<r<R„,  and  n°(r)=0,  otherwise.  Here  R, 
el)  l i.  e z 

is  the  outer  radius  of  the  E-layer  and  R^=R2~A  is  the  inner  radius. 

For  a thin  E-layer,  the  resulting  dispersion  relation  for  the  negative- 
mass  instability  can  be  expressed  as- 

, 2 
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(w-£u>  )‘ 
c 


/ fi.Aco  > 

2v  £ _ ( c 

yQ  l V 2c  / ’ 


(1) 


where  o>  is  the  complex  eigenf requency , is  the  electron 

cyclotron  frequency,  t is  the  azimuthal  harmonic  number,  and  use  has 
been  made  of  the  cold-fluid  limit  of  the  kinetic  dispersion  relation 
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for  the  negative-mass  instability.  In  Eq.  (1),  the  geometric  factor 
g is  defined  by  g=£^/(b_+b+) -k^RQ/(d_+d+) , and 
Ut(pR2)N^(pRc)-J^(pRc)Njl(pR2) 

V J1(FV  NI(P  V V (PR2  ) ’ 

R2PJi(pR2)Ni(pRc)-Jt(pRc)Ni(pR2) 


V " U£(pR2)¥Jl(pRc)-Ji(pRc)NJl(pR2) 


(2) 


b_-UJl(pR1)  / {RjPJ^pR^  ] , d_«l/b_, 

where  J (x)  and  N#(x)  are  Bessel  functions  of  the  first  and  second 
kind,  respectively,  the  prime  (')  denotes  (1/p)  (d/dr),  and  p is  defined 
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by  p = (w  /c  -k  ) , where  k Is  the  axial  wavenumber. 

The  nonlinear  consequences  of  the  negative-mass  instability 

can  be  investigated  within  the  framework  of  the  conservation 

equation  for  net  energy  balance. ^ We  Fourier  decompose 

perturbed  quantities  according  to  6A(^cft)3  £ 6A^  ^(r,t)exp(iJtO+ikz) , 

t,k  ’ 

where  k=  2irq/L  and  q is  an  integer.  Energy  conservation  can  then  be 


expressed  as 


t?  .V'  1;  ,irte«.k(r’t> 


where  k=|6gfc  k(r,  t)  1 2/&-+|6£e  k(r,t)|2/8n  is  the  perturbed  field 
energy  density,  and 

K - ~ [ dz [ de[  Cdrr [d3p(Y-l)mc2f  (x,p,t) 

L J-L/2  JO  J0  J 

is  the  average  electron  kinetic  energy  per  unit  length.  Here  L is  an 
effective  periodicity  length.  Within  the  context  of  a simple  quasilinear 
model,  the  spectral  energy  density  evolves  according  to^ 

!r€I,k<r-«)"2“5'S,k<r-t>  • <‘> 

i.  k 

where  the  growth  rate  = Imu)  is  determined  adiabatically  in  terms 

of  other  system  parameters  from  the  linear  dispersion  relation  in 
Eq.  (1).  Moreover,  for  a thin  E-layer  with  V^=*0  and  near-laminar 
flow  in  the  azimuthal  direction,  we  approximate  the  electron  kinetic 


energy  per  unit  length  by 


i-2n[  Cdrrn® 

Jo 


(r,t)  (y-l)tnc  , 


where  n^r.t)-/^^  !\ l/2  ls  the  avera8e  electron  density, 

and  ybc  »(l+r  w /cZ)X/  nc^  is  the  azimuthal  energy  associated  with  the 
c 

average  electron  motion  (V®“  ru^/y) . Here  u>c»eB0/mc  is  the  nonrelativistic 


r 
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electron  cyclotron  frequency. 

Substituting  Eqs.  (4)  and  (5)  into  Eq.  (3),  and  assuming  that 
the  electron  density  profile  maintains  a rectangular  shape  with 
variable  thickness  A(t),  we  obtain 

iNemc2F(A)^-  I (-5,kW£,k(0)exP<2^dfuJ*k(t')}  , (6) 

fR 

where  W . (0) h2tt I cdrre  . (r,t=0)  characterizes  the  energy  of  the  initial 
J Q 

field  perturbations,  Ne  is  the  number  of  electrons  per  unit  axial 
length,  and  F(A)>0  is  defined  by 


(2x.  [(l+x2)3/2-(l+x^)J/‘]  3x.(l+xf)A,M  £ 

f (*)-{— —~2— £3 XT if  * 

* ^2*^1  ^2  ' 

with  x,  a I R_/c  and  x,=w  (R0-A)/c.  Consistent  with  experimental 
z c 2 l c l 

observations,  in  obtaining  the  left-hand  side  of  Eq.  (6)  we  have 
approximated  dR2/dt^0  but  retained  terms  proportional  to  dA/dt. 


2*3/2, 


2*1/2 


(7) 


3.  NONLINEAR  EVOLUTION 


For  g>0,  it  is  evident  from  Eq.  (1)  that  the  linear  growth  rate 


t,k 


vanishes  once  the  layer  thickness  has  increased  to  the  critical 

value 


a 2c  ,2vg. 

V 


(8) 


Therefore,  making  use  of  Eq.  (6),  the  quasilinear  evolution  of  the 

system  can  be  summarized  as  follows.  At  t**0,  the  initial  thickness 

of  the  E-layer  is  Ai-A(t-0).  As  the  field  perturbations  grow,  A(t) 

continues  to  increase  {dA/dt>0  since  F(A)>0]  until  A(t-^»)“Aj  and  the 

instability  ceases.  (Here  we  assume  that  a single  1 value  is  excited.) 

Concomitant  with  the  increase  in  E-layer  thickness  is  an  increase  in 

8. 9 

azimuthal  energy  spread  of  the  electrons  composing  the  layer. 
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The  efficiency  n of  microwave  generation  is  defined  as  the  ratio 

of  electromagnetic  field  energy  generated  by  the  instability  to  the 

initial  kinetic  energy  of  the  E-layer.  In  the  ultrarelativistic 
2 2 

limit  with  Yq>>1,  X2»1,x^»1  and  A/Rq«1,  some  straightforward  algebra 
that  makes  use  of  Eqs.  (6)  and  (7)  shows  that  n can  be  approximated  by 


n 


(&f-A1)  . 


(9) 


For  future  reference,  we  also  introduce  the  average  power  per  unit 
axial  length  associated  with  the  microwave  production,  i.e., 

**i,k  * T jQ  rdr*E*.k(r,T)"Ei,k(r,0)  1 (10) 

where  T is  the  duration  of  the  radiation.  * 

Equations  (1)  and  (6)  have  been  solved  numerically  for  the  self- 

consistent  nonlinear  evolution  of  the  E-layer  thickness  A(t).  For 

present  purposes,  we  consider  transverse  electric  (TE  mode)  perturbations 
2 2 2 2 2 2 

with  k Rg=t  “a£n^Ro^Rc^  and  J£^a£,n^*^*  TyPical  results  are  summarized 

in  Figs.  2 and  3 for  Rc/Rg=1.5  and  Yq^A.S.  In  Fig.  2,  we  plot  normalized 

E-layer  thickness  A(t)/Rrt  versus  u>  t for  several  different  values  of 

u c 

2 2 

Budker’s  parameter  v=Nee  /me  . Evidently,  the  time-asymptotic 

thickness  of  the  E-layer  (A^)  is  a slowly  increasing  function  of  v 

(see  also  Eq.  (8)].  On  the  other  hand,  the  time  scale  required  for 

saturation  decreases  rapidly  with  increasing  v,  which  is  a consequence 

1 k 

of  the  fact  that  the  growth  rate  u^*  increases  with  v during  the 
initial  stages  of  instability  when  A<<Af  (Eq.  (1)J. 

We  note  from  Eqs.  (8)  and  (9)  that  the  efficiency  r\  of  microwave 
generation  is  also  an  increasing  function  of  Budker's  parameter  v. 

To  illustrate  this  point,  the  normalized  average  microwave  power  per 
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unit  length  P/w  (me  /e)  is  plotted  versus  v in  Fig.  3 for  R /R_»1.5, 

C c u 

Yq=4.5,  and  several  different  waveguide  modes  (l,n).  As  an  example, 

for  the  TE^  mode,  the  average  microwave  power  per  unit  length 

P increases  five-fold  when  v is  increased  from  0.07  to  0.15.  Also 

shown  in  Fig.  3 is  a plot  of  efficiency  n versus  v. 

4 . CONCLUSIONS 

Several  important  conclusions  follow  from  the  present  analysis. 

For  example,  the  time  evolution  of  the  E-layer  thickness  exhibits  a 

very  sensitive  dependence  on  Budker's  parameter  v [Fig.  2].  Moreover, 

the  average  power  and  efficiency  of  microwave  generation  are  enhanced 

considerably  by  an  increase  in  electron  density  [Fig.  3].  Introducing 

an  axial  energy  spread  into  the  present  analysis  reduces  the  growth 
1 k 5 9 

rate  substantially,  * thereby  slowing  the  rate  at  which  the 

E-layer  thickness  increases.  This  behavior  is  consistent  with  recent 
experimental  observations  by  Destler  et  al.  (cf. , Figs.  8 and  11  in 
Ref.  3). 
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FIGURE  CAPTIONS 

Fig.  1 Equilibrium  configuration  and  electron  density  profile. 

Fig.  2 Plot  of  A(t)/RQ  versus  u>ct  for  Rc/RQ=1.5,  y0=4.5,  A1/RQ- 

2xl0~2,  w . (0)/(N  roc2)=8xl0-3,  and  several  different  values 
of  Budker’s  parameter  v. 

- 2 2 

Fig.  3 Plot  of  normalized  microwave  power  per  unit  length  P/<Dc(mc  / e) 

versus  v [Eqs.  (10)  and  (6)]  for  R^/Rq^I.5,  Yq=4.5,  A^/Rq-2x10  , 

W (0)=8xl0-3  N me2,  and  several  different  modes  (i,n). 
ik  e 

Also  shown  is  a plot  of  efficiency  n versus  v. 


